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Abstract 

We present 0(a) YFS exponentiated results for wide angle Bhabha scattering 
at LEP/SLC energies using a new Monte Carlo event generator BHWIDE l.xx. 
Our calculations include two options for the pure weak corrections, as presented in 
Beenakker et al. and in Bohm et al. From comparison with the results of Beenakker 
et al. , Montagna et al. and Cacciari et al. , we conclude that the total precision 
of our BHWIDE results is 0.3%(0.5%) in the LEP1/SLC regime within ±100 MeV 
(+2.75/— 2.5 GeV) of the Z peak. For LEP2, the corresponding precision is currently 
estimated at 1.5%; the latter could be improved if the data in LEP2 so require. Both 
precision tags represent clear improvements over what is currently available in the 
literature. 
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1 Introduction 



As the final LEP1 data analysis begins and the initial stages of LEP2 materialize while 
the SLD prepares for the beginning of what may be its final phase, the subject of pre- 
cision calculations of wide angle Bhabha scattering becomes more and more interesting, 
particularly from the standpoint of realizing such calculations via a Monte Carlo event 
generator which would allow comparison between theory and experiment at the level of 
events in the presence of arbitrary detector cuts. Indeed, currently, the work of Beenakker 
et al. Jl| and of Montagna et al. represent the state of the art in the theoretical arena 
of the prespective calculations, and neither of these approaches has provided a genuine 
Monte Carlo event generator to allow comparison with arbitrarily cut experimental data. 
This has had the effect that the important Z physics parameter T eS has been extracted 
from the data using a subtraction of the contamination from the t-channel 7-exchange 
in the wide angle acceptance relevant to its measurement. A Monte Carlo event gen- 
erator with sufficient precision would obviate the need for this procedure, for example. 
Motivated by this and other such Z and LEP2 applications, we have developed an 0(a) 
Yennie-Frautschi-Suura (YFS) || exponentiated Monte Carlo event generator for wide 
angle Bhabha scattering at LEP1/SLC and LEP2 energies. We present this development 
in what follows. 

Specifically, our starting point will be the 0(a) YFS exponentiated Monte Carlo (MC) 
event generator BHLUMI l.xx developed by two of us (S.J. and B.F.L.W.) in ref. Q for low 
angle Bhabha scattering in the SLC/LEP luminosity regime 17 mrad < 9 e ^ < 150mrad, 
where 9 e ^ are the respective e, e CMS scattering angles. In order to arrive at an event 
generator valid for wide angles at both LEP1/SLC and at LEP2 to sufficient accuracy, we 
have had to introduce the effects of the Z exchange graphs into the calculations presented 
in ref. [|[] and we have had to introduce the effects of the pure weak one-loop corrections 
into those calculations as well. We stress that we have made contact with the pioneering 
0(a) MC program BABAMC of refs. |5], |6|] as well as with the semi-analytical program 
ALIBABA of ref. |TJ in that our pure weak corrections libraries are taken from these 
two cases: the user of our new wide angle Bhabha MC, BHWIDE l.xx, may choose 
which library he uses and we shall discuss both of them in what follows. Our exact hard 
bremsstrahlung amplitude, in the presence of the full set of s-channel and t-channel 7 and 
Z exchanges, we shall compute explicitly in what follows, using methods of the CALKUL- 
type |7J , as formulated by[| Xu et al. |J| . It is in this way that we have arrived at our new 
0(a) YFS exponentiated MC event generator BHWIDE l.xx which simulates realistic 
multiple photon radiative effects for wide angle Bhabha scattering in the LEP1/SLC and 
LEP2 energy regimes. 

We have compared our results with those of BABAMC, of the Monte Carlo integrator 
program SABSPV of ref. ||10||, of the semi-analytical program TOPAZO of refs. |TT|, 



of ALIBABA, of the MC event generator BHAGENE3 of refs. [13, 14 1, of the MC event 



generator BHAGEN95 of refs. [IIR pi TA [TJ, and of the MC event generator UNIBAB of 



ref. On this basis, we have arrived at a reliable estimate of the total precision of our 



We note that Kleiss and Stirling fel have introduced analogous methods as well. 
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results. The basis of this estimate is presented in detail in ref. [[HJ and we will review only 
its main features and the main features of the aforementioned attendant comparisons in 
what follows. 

More precisely, our BHLUMI l.xx Monte Carlo event generator realizes the process 
e+ipt) + e-(gi) — > e + {p 2 ) + e"^) + 71(^1) + . . . + ln{K) (1) 
via the YFS exponentiated cross section formula 



da 



^2aReB+2aB 



n=0 



3=1 J 



!)•••) ""n 



d 4 y 

d 3 p 2 d 3 q 2 
P2Q2 



(2) 



where the real infrared function B and the virtual infrared function B are given in refs. [§, 
^, 21, 22, 23| (for definiteness, we record them in the form we shall use presently, as this 
representation is not readily available in the current literature and is essential for practical 
applications of the type we pursue here), and where we note the usual connections 
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for the standard YFS infrared emission factor 
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(4) 



if Qf is the electric charge of / in units of the positron charge. Here, the ". . . " represent 
the remaining terms in S(k) obtained from the one given by respective substitutions of 
Qf, p\, Qf, q\ with corresponding values for the other pairs of the respective external 
charged legs according to the YFS prescription in ref. |§ (wherein due attention is taken 
to obtain the correct relative sign of each of the terms in S(k) according to this latter 
prescription) and in refs. ^5j, / ^ e, /' = / . We have explicitly the representations 



2a3t.B(p 1 ,q 1 ,p 2 ,q 2 ) + 2aB(p 1 ,q 1 ,p 2 ,q 2 ; k m ) = Rx{px,qx\ k r , 
+R 2 (pi,p 2 ; k m ) + R 2 (qi,q 2 ; k m ) - R 2 (pi,q 2 ; k r , 



Ri(p 2 ,q 2 ; K 
R 2 (qi,p 2 ] K 



with 
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1. 2 (A-5) 2 /A + uA . /A + u; 



where A = \f2pq + (p° — g ) 2 , u; = p° + g°, 5 = p° — q°, and fc m is a soft photon cut-off 
in the CMS (E^ < k m < £ feeam ). 

The YFS hard photon residuals in (EI), z = 0,1, are given in ref. [ffl for BH- 
LUMI l.xx so that this latter event generator calculates the YFS exponentiated exact 
0(a) cross section for e + e~ — > ee + n{^f) with multiple initial, initial-final, and final 
state radiation using a corresponding Monte Carlo realization of (0) in the low angle 
regime of the SLC/LEP luminosity monitor acceptance, where the e ± scattering angles 
8 e ± satisfy 17 mrad < 9 e ± < 150 mrad. In the next sections, we use explicit Feynman di- 
agrammatic methods and the results in refs. |TJ, ^ [J to develop the corresponding Monte 
Carlo realization of the respective application of (0) to the wide angle scattering regime, 
8 e ± > 150 mrad, for the Bhabha process e + e~ — > e + e~ + n(7). 

We turn first to the implementation of the required wide angle physics effects in 
the YFS hard photon residual /3q through 0(a). We follow this discussion with the 
corresponding analysis for the wide angle physics effects in the hard photon YFS residual 
f3\. In this way, we arrive at an exact 0(a) YFS exponentiated MC event generator valid 
for wide angle Bhabha scattering in the LEP1/SLC and LEP2 energy regimes. 

2 Born and O(a) contributions to (3q for wide angle 
Bhabha scattering at high energies 

In this section we develop the Born and O(o) contributions to the YFS hard photon 
residual /3q for the wide angle Bhabha scattering process. We use the low angle limit 
of the residual already presented in ref. | as a starting point and as a limiting case 
cross-check. 

Specifically, the hard photon residual /3q is defined as follows through 0(a): 



2' u dn dn ' 

where do x ~ loop '/ dQ is the differential cross section for wide angle Bhabha scattering com- 
puted through the 1-loop correction and da Born /dO, is the respective Born differential 
cross section. In ref. [HI, the right-hand side of (|8]) was computed in the low angle regime. 
Here, we need to compute (§) in the wide angle scattering regime. We do this as follows. 
First, we note that the exact expression for the wide angle scattering Born differential 
cross section in (P) is well-known 
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where | _Axt 1 2 is a spin averaged lowest order matrix element squared, given e.g. in eq. (22). 
Secondly, the complete 1-loop corrected wide angle Bhabha scattering differential cross 
section in (H) is also known and we use two different versions of it in our work, one taken 
from ref. [[J and one taken from ref. where it is generally accepted that the latter 
version, which is the more recent of the two, is in fact the more up-to-date of the two [pOfl 
(and hence it is the option which we illustrate in our comparisons in Sect. 4). In our 
Monte Carlo event generator these two versions for da 1 ^ loop / dfl correspond to two choices 
for our electroweak library module, one of which the user specifies in his input file ||26|1 . 
It is this way that we have realized the wide angle Bhabha scattering YFS hard photon 
residual (3q through 0(a). 

We turn next to the wide angle Bhabha scattering YFS hard photon residual (3\. This 
we do in the following section. 

3 The YFS hard photon residual f3\ for wide angle 
Bhabha scattering at high energies 

In this section we present our determination of the 0(a) YFS hard photon residual j3\ 
for wide angle Bhabha scattering at high energies. We start with the defining equation 
for \. 

The hard photon residual Pi, to 0(a), is defined as follows 

where da B1 / kdkdQ^dQ is the respective O(a) bremsstrahlung differential cross section 
into the solid angles dVt 1 and dQ for the photon and positron respectively when the 
photon energy lies between k and k + dk. Thus, to specify fix we need to specify the hard 
bremsstrahlung differential cross section on the RHS of (|10|). We now turn to this. 

Specifically, the hard bremsstrahlung differential cross section required in fllOl ) can be 
obtained via the standard methods from the corresponding helicity amplitudes M. ( A e + , A e - , 
A' e+ , A' e _, A 7 ), where Ay , / = e + , e~ is the incoming (outgoing) fermion helicity and A 7 is 
the photon helicity. We have 

dcr B1 1 (n ) 2 



kdkdVt^dVt 512vr 5 s 2E b (E b - k) 
with 



\ M \ 2 = \ E IMA e+ ,A e -,A; + ,A^,A 7 )| 2 (12) 



e i7>j'= 



being the spin averaged squared matrix element. The formula ([TTJ) is given in the CMS 
of the incoming beams, where E b is the beam energy. Thus, our determination of f3\ will 
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be completely specified when we give our results for A4(X e +, A e -, X' e+ , A' e _, A 7 ). We now 
turn to this. 

More precisely, using the methods of ref. ||, we get the following representation for 
the non-vanishing helicity amplitudes: 
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where our notation is as follows 
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Invariants: 



s = 2p 1 q 1 , s' = 2p 2 q 2 , 
t p = -2p ± p 2 , t q = -2gig2, 

u = -2piq 2 , u' = -2q x p 2 . (14) 



Electroweak couplings: 



G G 
a L = — (v e + a e ), a R = — (v e -a e ), (15) 



where 



G= . 7i ' a e = , ■ a ^— , w e = a e (l -4sin 2 ^). (16) 

sin Vw cos &w 4 sin Uw cos Uw 

Spinor functions: 

< pq > 



F(p, q) = 



< pq >* 



G fo^*) = <r < k><ks> > (17) 

where p,q,r,s denote fermion four-momenta in the massless limit, while k denotes 
the photon four-momentum, and the * stands for complex conjugation. Spinor 
products < .. > are defined as 

< pq >= Vp^e^ - y/p^qle*', (18) 
where p± = p ± p z , p± = p x + ip y = ^/p+P- e 1 ^. 
Propagator factor: 

R ^ = l + x -MlZ( X )M z V z < 19 » 



Specific notation: 



F 1 = F(p 1 ,q 1 ), F 2 = F(p 2 ,q 2 ), F 3 = F(p 1 ,p 2 ), 

F 4 = F( qi ,q 2 ), F 5 = F( Pl ,q 2 ), F 6 = F(q 1 ,p 2 ), 

Gi = G(p u p 2 ,q u q 2 ;k), G 2 = G(q 2 ,q 1 ,p 1 ,p 2 ;k), 

G 3 = G(q 2 ,p 2 ,p l ,q l ;k), G A = G(p 1 ,q 1 ,p 2 ,q 2 ;k), 

c = i2V2e 3 . (20) 
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The above helicity amplitudes have been obtained in the massless fermion approxima- 
tion. However, to get a precise description of the photon radiation over the whole phase 
space (particularly for collinear configurations) the finite fermion masses have to be taken 
into account. This can be accomplished by adding to the matrix element \Ai\ 2 of eq. flT2]) 
the mass correction term [f7| 



,2 ™2 

|2 



5\M\ 2 mc = -e< 



L(^ i ^ o(s/ '^' m/)|2 + (^ i>(,,(s ' / ''-' /)|: 



m l \ \ A { „ 4- . M2 , m l \KA + „./M2 



(/cp 2 ) 2 ' ' (kq. 



2 



(21) 



where 



|M,(«> t, «)| 2 = e 4 |^ [|i?( s , 4)| V + |i2(s, a 2 ,)] V + 2\R{s, a L a R )\H 2 ] 
+^ [|i2(t, 4) I V + |i?(t, a|)|V + 2\R(t, a L a R ) \ 2 s 2 ] 
+^2U[R*(s,al)R(t,al) + R*(s,a 2 R )R(t,a 2 R )]^ (22) 

is the lowest order matrix element. 

We have also used the results of ref. |7j for the matrix element \A4\ 2 as a cross check 
and we have found that the two sets of results are in very good agreement with one 
another, well below the desired technical precision of 0.01% of our analysis for example. 
We also note that an equivalent representation of the results in ref. [|7| has been given by 
Kleiss in ref. f27j . On introducing these two sets of results into the formula (|T0"D for f3\ 
and implementing the resulting expression into our YFS Monte Carlo program for Bhabha 
scattering, we arrive at the Monte Carlo event generator BHWIDE 1.00 for wide angle 
Bhabha scattering at high energies^. 

We will now illustrate the application of BHWIDE 1.00 to LEP1/SLC and LEP2 
physics scenarios in the next section. 



4 Results and comparisons 

In this section, we present sample Monte Carlo data which we use to compare our pre- 
dictions from BHWIDE 1.00 to those of related calculations as reported in ref. | 20| . We 
discuss both LEP1/SLC energies and LEP2 energies. 

Specifically, for our comparisons we use the same event selection (ES) cuts as defined 
in ref. |p0[| : 



2 We need to stress that, in implementing the results (|lT|) — (22) and the corresponding results from 
ref. into (|l0[), care must be taken to compute the two terms on the RHS of ([ll]) in the same manner 
with the same massive (massless) limits of the corresponding parts of each respective term so that the 
result for /3i is numerically stable; we have done this. 
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No. 




TOPAZ0 


BHWIDE 


BHAGENE3 


ALIBABA 


BHAGEN95 










(a) BARE acol max 


= 10° 










1. 


88.45 


0.4579 ± 


.0003 


0.4560 ± .0004 


0.4495 ± 


0016 


0.4575 ± 


.0003 


0.4578 ± 


.0002 


2. 


89.45 


0.6452 ± .0002 


0.6429 ± .0006 


0.6334 ± .0023 


0.6440 ± .0003 


0.6445 ± .0003 


3. 


90.20 


0.9115 ± 


.0002 


0.9087 ± .0008 


0.8997± 


0033 


0.9090 ± 


.0004 


0.9095 ± 


.0004 


4. 


91.19 


1.1846 ± 


.0002 


1.1797 ± .0010 


1.1847± 


0033 


1.1840 ± 


.0004 


1.1822 ± 


.0005 


5. 


91.30 


1.1639 ± 


.0002 


1.1592 ± .0009 


1.1667 ± 


0033 


1.1636 ± 


.0005 


1.1619 ± 


.0005 


6. 


91.95 


0.8738 ± 


.0002 


0.8711 ±.0007 


0.8856 ± 


0028 


0.8769 ± 


.0003 


0.8742 ± 


.0004 


7. 


93.00 


0.4771 ± 


.0002 


0.4761 ±.0005 


0.4808 ± 


0019 


0.4814 ± 


.0001 


0.4796 ± 


.0002 


8. 


93.70 


0.3521 ± 


.0002 


0.3512 ± .0004 


0.3521 ± 


0013 


0.3556 ± 


.0001 


0.3550 ± 


.0001 










(b) BARE acol max 


= 25 u 










1. 


88.45 


0.4854 ± 


.0003 


0.4808 ± .0005 


0.4699 ± 


0016 


0.4833 ± 


.0003 


0.4826 ± 


.0002 


2. 


89.45 


0.6746 ± 


.0003 


0.6699 ± .0006 


0.6593 ± 


0023 


0.6727 ± 


.0003 


0.6710 ± 


.0003 


3. 


90.20 


0.9438 ± 


.0003 


0.9387 ± .0008 


0.9279 ± 


0033 


0.9425 ± 


.0003 


0.9384 ± 


.0004 


4. 


91.19 


1.2198 ± 


.0003 


1.2130 ± .0010 


1.2169 ± 


0034 


1.2187± 


.0004 


1.2133 ± 


.0005 


5. 


91.30 


1.1989 ± 


.0003 


1.1924 ± .0010 


1.1995 ± 


0034 


1.1982 ± 


.0004 


1.1928 ± 


.0005 


6. 


91.95 


0.9054 ± 


.0002 


0.9011 ±.0007 


0.9124± .0026 


0.9089 ± 


.0003 


0.9014 ± .0003 


7. 


93.00 


0.5040 ± .0002 


0.5013 ± .0005 


0.4996 ± .0019 


0.5054 ± .0002 


0.5027 ± .0002 


8. 


93.70 


0.3777 ± 


.0002 


0.3749 ± .0004 


0.3689 ± 


0013 


0.3782 ± 


.0001 


0.3771 ± 


.0001 



BARE acoL 



10° 



BARE acoL 



25° 



0.02 
0.01 
0.00 
-0.01 
-0.02 



Q--Q-REF 
"■REP 



❖ BHWIDE 

* BHAGENE3 
O ALIBABA 
□ BHAGEN95 
X TOPAZO a R ~L 



-x- 



,_ r . 



® 
□ 



No. of Energy point 
i i i 



0.02 
0.01 
0.00 
-0.01 
-0.02 



o 



r-Q- REP 
"■REP 



BHWIDE 
★ BHAGENE3 
O ALIBABA 
□ BHAGEN95 
X TOPAZO 



<> 



$ 



. ® . 



-x- 



-x- 

...a... 

<l> 



„No. of Energy point 
i i , , , , 



1. 2. 3. 4. 5. 6. 7. 8. 1. 

Figure 1: Monte Carlo results for BARE trigger, for two values (10° and 25°) of acollinearity cut. 
Center of mass energies (in GeV) close to Z peak. In the plots cross section cjref from TOPAZO 
is used as a reference cross section. Cross sections in nb. Two horizontal dotted lines indicate the 
0.3% band, for reference. 



• BARE ES - we require 40° < e - < 180°, 
of 10°, 25° and E min = 1 GeV for both e~ 

• CALO ES - we require the same cuts as in 



0° < 6 e + < 180°, with acollinearity cuts 
and e + ; 

the BARE case but with E min = 20 GeV 
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No. 




TOPAZ0 


BHWIDE 


BHAGENE3 


UNIBAB 


BHAGEN95 










(a) CALO acol max 


= 10° 










1. 


88.45 


0.4533 ± 


.0004 


0.4523 ± .0004 


0.4467 ± 


0008 


0.4490 ± 


.0010 


0.4524 ± 


.0001 


2. 


89.45 


0.6387 ± .0004 


0.6377 ± .0006 


0.6302± .0011 


0.6358 ± .0012 


0.6370 ± .0002 


3. 


90.20 


0.9023 ± 


.0003 


0.9016 ± .0008 


0.8920 ± 


0015 


0.9021 ± 


.0014 


0.8990 ± 


.0003 


4. 


91.19 


1.1725 ± 


.0001 


1.1707 ± .0010 


1.1767± 


0021 


1.1772 ± 


.0016 


1.1689 ± 


.0004 


5. 


91.30 


1.1520 ± 


.0001 


1.1505 ± .0009 


1.1571 ± 


0020 


1.1559 ± 


.0016 


1.1491 ± 


.0004 


6. 


91.95 


0.8649 ± 


.0001 


0.8646 ± .0007 


0.8795 ± 


0015 


0.8689 ± 


.0012 


0.8660 ± 


.0003 


7. 


93.00 


0.4723 ± 


.0001 


0.4725 ± .0005 


0.4796 ± 


0008 


0.4733 ± 


.0008 


0.4761 ± 


.0001 


8. 


93.70 


0.3486 ± 


.0001 


0.3486 ± .0004 


0.3507± 


0006 


0.3486 ± 


.0007 


0.3526 ± 


.0001 










(b) CALO acolmax 


= 25° 










1. 


88.45 


0.4769 ± 


.0004 


0.4742 ± .0004 


0.4696 ± 


0008 


0.4733 ± 


.0010 


0.4741 ± 


.0001 


2. 


89.45 


0.6638 ± 


.0003 


0.6615 ± .0006 


0.6556 ± 


0011 


0.6619 ± 


.0012 


0.6599 ± 


.0002 


3. 


90.20 


0.9297 ± 


.0003 


0.9278 ± .0008 


0.9207± 


0012 


0.9302 ± 


.0014 


0.9236 ± 


.0003 


4. 


91.19 


1.2025 ± 


.0003 


1.1994 ± .0010 


1.2074± 


0021 


1.2073 ± 


.0016 


1.1954 ± 


.0004 


5. 


91.30 


1.1819 ± 


.0003 


1.1790 ± .0010 


1.1879 ± 


0021 


1.1860 ± 


.0016 


1.1753 ± 


.0004 


6. 


91.95 


0.8924 ± 


.0003 


0.8909 ± .0007 


0.9058 ± .0016 


0.8965 ± 


.0012 
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Figure 2: Monte Carlo results for CALO trigger, for two values (10° and 25°) of acollinearity cut. 
Center of mass energies (in GeV) close to Z peak. In the plots cross section cjref from TOPAZO 
is used as a reference cross section. Cross sections in nb. Two horizontal dotted lines indicate the 
0.3% band, for reference. 



for the final 'fermion' energy which is the e~(e + ) energy if the there are no photons 
nearby and which is the e~(e + ) plus the photon energy if the photon is within a 
cone of half-angle 1° from the e~(e + ), respectively. 
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24.976 ± .115 


35.740 ± .080 
30.270 ± .070 
25.960 ± .060 
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30.296 ± .016 
25.958 ± .014 




(b) CALO acol max = 25° 
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3. 


39.741 ± .049 
33.698 ± .042 
28.929 ± .036 


40.487 ± .025 
34.336 ± .017 
29.460 ± .013 


39.170 ± .280 
32.400 ± .190 
27.840 ± .160 


39.521 ± .158 
33.512 ± .135 
28.710 ± .116 


40.240 ± .100 
34.100 ± .080 
29.280 ± .070 
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34.287 ± .018 
29.409 ± .015 



CALO acoL 



10° 



CALO acoL 



25° 



0.10 



0.05 



0.00 



-0.05 



-0.10 



■I i i r 

O TOPAZO 



<r — ^ref ^ BHAGENE3 
o- REF O UNIBAB 

□ BHAGEN95 

• SABSPV 

X BHWIDE C„ 



-4- 



No. of Energy point 
i i i i. 



1. 



2. 



3. 



0.10 



0.05 



0.00 



-0.05 



-0.10 



T 



1 1 1 1 1' 

O TOPAZO 



c7— £7rep w BHAGENE3 
<r REF O UNIBAB 

□ BHAGEN95 

• SABSPV 

X BHWIDE C„ 



+ 



No. of Energy point 
i i i i. 



1. 



2. 



Figure 3: Monte Carlo results for CALO trigger, for two values (10° and 25°) of acollinearity cut. 
Center of mass energies close to VF-pair production threshold {Ecu'- 1- 175 GeV, 2. 190 GeV, 
3. 205 GeV). In the plots cross section oref f r ° m BHWIDE is used as a reference cross section. 
Cross sections in pb. Two horizontal dotted lines indicate the 1.5% band, for reference. 



For the BARE ES acceptance cuts, we show in fig. || the comparison of the BHWIDE 
results with those of ALIBABA, BHAGEN95, BHAGENE3 and TOPAZO, where for def- 
initeness we plot the ratio {a a — &ref) / &ref for each calculation A, A = ALIBABA, 
BHAGEN950, BHAGENE3, BHWIDE, and TOPAZO, using the TOPAZO cross section as 
<tref- This we do for the 8 CMS energy values: 88.45, 89.45, 90.20, 91.19, 91.30, 91.95, 
93.00, 93.70 GeV, which are denoted in the figure as energy points 1,2,... ,8, respec- 
tively. We see that BHWIDE agrees with the semi-analytical programs, ALIBABA and 
TOPAZO, to within a few per mille at the Z peak whereas off the peak BHWIDE remains 
generally within 1% of the semi- analytical programs. The difference between BHWIDE 
and the Monte Carlo program BHAGENE3 is also at the few per mille level at the Z peak 

3 Here we use recently updated results of BHAGEN95, obtained from the authors. 
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but it is as much as 2.27% off the peak. The agreement between BHWIDE and BHA- 
GEN95 is also at the few per mille level at the Z peak; off the peak, two programs remain 
within 1% of one another, with the better agreement holding for the looser acollinearity 
cut of 25°. These results give us a handle on the total precision of BHWIDE in the Z 
resonance regime, as we discuss presently. 

Continuing in this way, we show in fig. ^] the similar type of comparison of the predic- 
tions of the programs with the CALO ES. In fig. ALIBABA is no longer applicable (it 
does not handle the CALO ES) and UNIBAB appears — it was too slow to participate 
in the BARE ES comparisons. Again, we use TOPAZO for the reference cross section and 
we plot the same ratio [a a — &ref) / &ref for the same 8 energy points as in fig. |l[ At the 
peak, BHWIDE is within a few per mille of TOPAZO and BHAGEN95 and it is within 7 
per mille of UNIBAB and BHAGENE3. Off peak, BHWIDE remains within 6 per mille 
of TOPAZO whereas it remains within 0.8% of UNIBAB, within 1.2% of BHAGEN95 and 
within 1.8% of BHAGENE3; the better agreement with the latter three programs holds 
for the looser acollinearity cut whereas for TOPAZO the situation is reversed. Based on 
these and related comparisons as described in ref . ^(J , including both the results and the 
physics approximations in the various programs as presented in this last reference, we 
conclude that, for the CALO ES, within ±100 MeV of the Z peak, the total precision 
of BHWIDE is 3 per mille and off peak, within +2.5/ — 2.75 GeV thereof, we set this 
precision at 5 per mille in the LEP1 energy regime. For reference the three per mille 
band is indicated by the two horizontal dotted lines in figs. [I] and |2|. This precision tag 
should be compared to that for the only other published wide angle Bhabha scattering 
Monte Carlo event generator in wide use at LEP/SLC, namely BABAMC ||, whose total 
precision on pure QED was set at 1% in the Z peak region in refs. |2"5fl . (In practice, 
ALIBABA and/or TOPAZO would be used to determine the non-QED and higher order 
QED part of the respective cross section.) 

Turning next to LEP2 energies, we show in fig. |3| the comparison of the results of the 
six programs BHAGEN95, BHAGENE3, BHWIDE, SABSPV, TOPAZO, and UNIBAB 
in the same format as in figs. [I], 0, using in this figure BHWIDE for the reference cross 
section (Jref for the three CMS energy points 175, 190, and 205 GeV . We use the CALO 
ES only here. SASPV and BHAGEN95 are within 2% of BHWIDE for both acollinearity 
cuts, with SASPV within 1.5% of BHWIDE and with BHAGEN95 maximally at 0.55% 
above SASPV. TOPAZO and UNIBAB deviate from BHWIDE by as much as 1.9% and 
2.4% respectively, with the worse agreement holding for the looser (tighter) acollinearity 
cut, respectively. These deviations are consistent with the leading logarithmic accuracy 
one expects for these two Z peak optimized codes. Similarly, the deviations between 
BHWIDE and BHAGENE3 by as much as 3.8% as shown in fig. |3| are also consistent with 
what one expects from the leading logarithmic Z peak optimized nature of BHAGENE3. 
On the basis of these and related comparisons, we estimate the total precision of BHWIDE 
at the LEP2 energies as 1.5%, conservatively. For reference, the 1.5% band is indicated 
by the two horizontal dotted lines in fig. ^. 

We end this section by noting that, in addition to the comparisons just presented, we 
have also checked that the pure QED 0(a) predictions of BHWIDE are within 0.05% of 
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those of the exact 0(a) MC OLDBIS of Ref. \22\ - a modernized version of the MC 
of Ref. []3(J . This gives us additional confidence in the technical component of our total 
precision estimate for BHWIDE 1.00. 



5 Conclusions 

In this paper, we have presented a new multiple photon Monte Carlo for wide angle 
Bhabha scattering at LEP1/SLC and LEP2 energies in which the respective multiple 
photon effects are realized on an event-by-event basis and in which the infrared singulari- 
ties are cancelled to all orders in a via YFS exponentiation. This Monte Carlo calculation 
contains the exact 0(a) result and features two choices for the respective pure weak cor- 
rections at 0(a), that in ref. [Q and that in ref. 0. It thus corresponds to the exact O(o) 
YFS exponentiated treatment of wide angle Bhabha scattering. 

We have illustrated our new calculation at both LEP1/SLC energies and at LEP2 ener- 
gies, for two types of event selection, the BARE and CALO selections of ref. [2D|, which fea- 



ture two choices of the acollinearity cut, 10° and 25°. In our illustrations, we compared our 
predictions with those of refs. [p], pi], 0, [y], [18|, [19| . We found in general a good agreement 



of the various calculations, an agreement consistent with the levels of approximations and 
realms of applicability of the respective codes. In this way, we arrived at the precision tags 
of 0.3% for BHWIDE at the Z peak and of 1.5% at LEP2 energies. The program is avail- 
able from the authors at the WWW URL [http:/ /enigma.phys.utk.edu/pub/BHWIDE 



In summary, we have developed a new Monte Carlo event generator for wide angle 
Bhabha scattering at LEP1/SLC and LEP2 energies in which the infrared singularities 
are cancelled to all orders in a via YFS exponentiation of the respective multiple photon 
effects. We look forward with excitement to its application to LEP1/SLC and LEP2 data 
analyses. 
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